Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



Implications of the HST/FGS parallax of SS Cygni 
on the disc instability model 



M.R. Schreiber^'^ and B.T. Gansicke^ 
Universitats-Sternwarte, Geismarlandstr.il, D-37083 Gottingen, Germany 

Department of Astronomy, University of Cape Town, Private Bag, Rondebosch 7700, Cape Town, South Africa 



Received 



accepted 



Abstract We analyse the consequences of the recently measured parallax of SS Cygni ( Harrison et al. 1999| ) on 
the accretion disc limit cycle model. Using the observed long term light curve of SS Cyg and d = 166 ± 12pc, we 
obtain for the mean mass transfer rate Adu — 4.2± 1.7 x lO^^gs"^. In addition, we calculate the vertical structure 
of the accretion disc taking into account heating of the outer disc by the stream impact. Comparing the mean 
accretion rate derived from the observations with the calculated critical mass transfer rate, we find that the disc 
instability model disagrees with the observed long term light curve of SS Cyg as Mtr is greater or similar to the 
critical mass transfer rate. The failure of the model indicated by this result can be confirmed by considering that 
the accretion rate at the onset of the dechne should be exactly equal to the value critical for stability. In contrast 
to this prediction of the model, we find that the accretion rate required to explain the observed visual magnitude 
at the onset of the dechne must be significantly higher than the critical mass transfer rate. Our results strongly 
suggest that either the usually assumed temperature dependence of the viscosity parameter a is not a realistic 
description of the disc viscosity, that the mass transfer rate in SS Cyg noticeably increases during the outbursts 
or, finally, that the HST distance of 166 ± 12 pc, is too high. 
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1. Introduction 



In non- magnetic cataclysmic variables (CVs) (Warner 



|1995 , for an encyclopaedic review) a white dwarf accretes 
material from the Roche-lobe filling secondary via an ac- 
cretion disc. In many CVs, sudden brightenings - called 
dwarf nova outbursts - are observed. These eruptions are 
generally thought to result from thermal instabilities in 
the accretion discs (the disc instability model (DIM), see 
e.g. p4uiil^^o 1993b| ; pudwlg eL al. 1994| , for parameLcr 



the cold state, ahot and acoid, with ahot/ctcoid typically 
5-10 but with unexplained exceptions ahot/o^coid ~ 10^ 
e.g. for WZ Sge. Although the change of a is a some- 
what arbitrary assumption, the DIM has become the gen- 
erally accepted explanation of dwarf nova outbursts as 
it a grees with obser vational and theoretical constraints 
(see Cannizzo 1993b ). Recently, additional ef fects such as 
disc radius variations ( Hameury et al. 1998 ) , irradiation 
( Hameury et al. 199E ; Schreiber et al. 2000a|; ^chreiber fc 



vertical dlyc yti'iicLuro for a glveu radius show a hysLere 
sis relaiLloii 111 Llie stirface-deiislLy vy. accreLloii rale plaue 



Gansicke 2001), mass transfer v ariations ( Schreiber et al 



the iiu-called "S-curve'\ lu prluelple, a UmU cycle be- — in the model 



studies ). CalcmaLloiis of Llie yLeady-sLale yomLlou of the ^qqq^^ [suat-Menard et al. 200l|) and the effects of the ac- 



cretion stream impact (Schreiber & Hessman 1998; Buat 



Menard et al. 2001 ; gtehle et al. 2001 ) have been included 



haviour can be maintained by a mean mass-transfer rate 
which is between the critical rate needed to stay in a cold 
quiescent state and that needed to remain in the hot out- 
burst state. However, the obtained light curves using a 
constant viscosity parameter a show only low amplitude 



outbursts (Smak 1984). In order to match the observed 
outburst light curves, it is necessary to assume a drastic 
change of the viscosity parameter a between the hot and 



Send offprint requests to: M.R. Schreiber, email: mschrei@uni- 
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2. The absolute brightest dwarf nova: SS Cyg 

The dwarf nova SS Cyg is among the visually brightest 
CVs and has a detailed long term light curve covering 
more than a century. SS Cyg has, therefore, become the 
reference object for accretion disc instability analyses, and 
its observed light cu rve has been used to calibr ate the vis- 
cosity parameter a. Cannizzo fc Mattel (1992 ) examined 
the long term light curve of SSCyg (Fig. 1) in detail in 
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j Figure 2. The absolute magnitudes of dwarf novae 
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Figure 1. The visual long term light curve of SSCyg 
compiled from observations made by members of 
the Association Francaise des Observateurs d'Etoiles 
Variables. 



order to test the limit cycle model (Cannizzo 1993a) and 
found a mean outburst duration of tout — 10.76 d and a 



adapted from Warner (1987) corrected to a mean incli- 
nation of i = 56°. 7 and with mv — 8.5 for the maxi- 
mum magnitude of SS Cyg. Due to the increased distance, 
SS Cyg has become the absolute brightest dwarf nova. 



and its maximum absolute magnitude is brighter than one 
would expect from its orbital period ( Cannizzo 199^ ). 

The mean accretion rate in SS Cyg inferred from its 
long term light curve depends, in addition to the absolute 
magnitude, on the inclination i of the accretion disc. This 
system parameter is somewhat uncertain, most studies of 



mean cluiesceiice Lime of iqui = 38.71 d (coiTeHpondliig to 
a meairnycle dut'atloii of tcyc = 49.47 d). In a latet' paper. 



SS Cyg suggest i = 37° ± 5° ( |Ritter fc Kolb 1998|). A sig- 
nificantly higher value of z ~ 57° was derived by Voloshina 



Khruzina (1997 ), but, throughout this paper we assume 



the lower value. To account for the influence of the incli- 



[Canniz^zo fc Mattel ( 1998p again examined the light^ curve ^^^-^^ absolute magnitude we follow [Paczynski fc 



of SS(j^yg, derlvliig Llie rl^e aiid decay Llme^. Oiily the gchwarzenberg-Czerny (1980|): 



mean values are Important in the context of oui' analysis: 
iris = 0.5 d, tdoc — 2.5 d. These values are obtained using 
a conservatively low estimate of the outburst magnitude 
of m-u = 8.5. 

However, even though the brightness variation of 
SS Cyg has been monitored to a high level of preci- 
sion, its distance has been highly uncertain. As a conse- 
quence, the absolute magnitude and accretion rates, have 
also been ill-defined. Recently, a high-quality parallax of 
SS Cy g was obtained using the Hubble Space Telescope 



AMv(i) = -2.5 log 



(0.4 -t- 0.6 cos i) cos i 
OA 



(1) 



(HST) [Fine Guidance Sensor (FGS), which moves the sys- 
tem to a distance 1.5 — 2 times larger than previously 
thought ([Harrison et al. 1999|). The HST/FGS distance 



where we used a limb darkening coefficient u = 0.6. Notice, 
Eq. (1) gives the inclination correction corresponding to a 
mean inclination i — 56°. 7 (Fig. 2). 

Throughout this paper we use the system parameters 
of SS Cyg given in the literature, i.e. the primary mass 
_Mwd = 1.19±O.O5M0 ( iRitter fc Kolb 1998| ; Teescling 
1997), secondary mass Algcc — 0.70 AIq, and orbital pe- 
riod in hours Phr = 6.6 ( [Ritter fc Kolb 1998| ). 



of 166 ± 12 pc is confirmed by Beuermann (2000), who 
showed that the accretion disc in SS Cyg contributes sig- 
nificantly to the observed K band flux. As a consequence 
of the now firmly established (larger) distance, the ab- 
solute magnitude of SS Cyg in outburst increased, too. 
Figure 2 shows the absolute visual magnitudes of dwarf 
novae at outburst maximum as a function of their orbital 



period, adapted from Warner (1987 ). The arrow displays 
the impact of the HST/FGS parallax on SS Cyg: the sys- 
tem is now the absolute brightest dwarf nova in outburst 



3. Deriving the critical and the mean mass 
transfer rate 

The goal of this section is to check if the "standard" disc 
instability model can reproduce the absolute magnitude 
of the outbursts of SS Cyg, taking into account the larger 
HST/FGS distance, the inclination, and the mean out- 
burst time scales derived from the long term light curve. 
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Table 1. Assumed system parameter 



Binary parameter Outburst parameter 





= 1.19 ±0.05 Mq 


^out 


10.76 d 


Msec 


0.70 Mq 


iqui — 


38.71 d 


i 


37° ± 5° 


icyc ~ 


49.47 d 


P 


= 6.6 hr 


^ris 


0.5d 


d 


166±12pc 


^dcc — 


2.5d 




= 0.90 ± 0.05i?Li 


mv = 


8.5 ±0.1 



Stream impact and irradiation heating efficiencies 
77, = 1.0 ±0.25 (1 - /3) = 0.25 ± 0.25 



Contribution of the irradiated secondary 

'TlV,scc > 11 



Figure 3. Simulated light curves using the orbital param- 
eters of SSCyg, ah = 0.5, = 0.05, log(Mtr)y = 17.7 
(upper panel) and log (Mtr) = 17.4. 

3.1. The mean mass transfer rate 

Figures shows a light curve computed with the FE/FD 



layer using Eq. (4) and (5) of Schreiber et al. (2000t) with 
an albedo of the accretion disc of /3 = 0.75 ± 0.25. To 
account for the contribution of the irradiated secondary 
(K4/5) we used the code BINARY+ ± written by A. van 
Teeseling (van Teeseling et al. 1998) and found even for 



an extremely high accretion rate of Mout = 4 x lO^^gs" 
the visual magnitude of the irradiated secondary being 
™v,scc ^ 11- We consider sec = 11 as an upper limit 
for the contribution of the irradiated secondary through- 
out this paper. 

The accretion rate necessary to reproduce the observed 
visual flux during the outburst can be described with a 
power-law: 



DIM code used for the first time in Schreiber & Hessman 
(1998|)]and described in detail in [Schreiber et al. (2000b|) 



M„„, = 3.62 X 10^^ g s 



\ 166pc 



We adapted the parameters suggested by Ludwig et al 
(1994|)]and two different values for the mass transfer rate. 

During an outburst cycle the disc subsequently goes 
through a phase of quiescence, rise, adjustment, quasi- 
stationary outburst state, and decay. After the heat- 
ing front has brought the entire disc into the hot state 
(rise), the disc redistributes its mass (adjustment) until 
it accretes without radial variations of the accretion rate 
(quasi-stationary state). Almost all the mass accreted dur- 
ing an outburst cycle is transferred onto the white dwarf 
during this quasi-stationary outburst state (Fig. 3). Thus, 
to derive a lower limit for the mass accreted during a 
mean outburst cycle, one simply has to multiply the ac- 
cretion rate of a stationary standard disc which reaches 
mv = 8.5 ± 0.1 with the mean quasi-steady state dura- 
tion, ~ 5.76 dQ. Throughout this paper we assume a ^out = (3.62 ± 1.47) x 10'* gs 
very large disc for SS Cyg during outburst, i.e. (90 ± 5)% 



/ J \ 3.990 / . X 1.499 



37= 



, r \ -1.449 , 



1.0 



1.19 M( 



'J0 



-0.023 /i^^Q 



5.83 



-1.130 



0.058 



0.75/ V8.5/ 



-15.508 



(3) 



With the parameters and their uncertainties given in 
Table 1 . , this results in 



18 „o-l 



(4) 



of the primary's Roche radius i?Li- (Ritter 1980; Barrop- 
Allin ^ Warner 19"96| ) (0.9 Rli = 5.83 x 10'" cm for the bi- 



The accreted mass has to be transferred to the disc during 
a mean outburst cycle by a mean mass transfer rate: 



nary parameter of SSCyg). Following Buat-Mcnard et al 
(200l|)] we describe the additional heating of the outer disc 
due to the stream disc impact with 



M^^ ^ i^Wqs ^ ^4 2± 1.7) X IQi^gs-i. 

^CVC 



(5) 



GMwdAftr 1 



i?d ^ R 



(2) 



= '^^ 2i?d 27ri?dAi?h/''P V^K 

where i?d is the outer disc radius and rji the efficiency of 
disc heating. In this section we assume 771 = 1.0± 0.25 and 
the Keplerian kinetic energy being released over a radial 
extent of A i?hs = 0.1 i?d with an exponential attenuation 
(Buat-Menard et al. 2001). In addition, we include irra- 
diation of the disc by the white dwarf and the boundary 



This value has to be compared with the critical mass trans- 
fer rate. 

3.2. The critical mass transfer rate 

In the absence of stream impact heating the critical accre- 
tion rate is given by the accretion rate at the upper turning 
point of the S-curve Ma and can be approximated by a 
power-law, i.e. 



Assuming an adjustment time of ta 



2d. 



Ma = Af< 



crit 



9.5 X lO^^gs^^i?: 



,2.64 
10 



wd 



-0.88 



(6) 
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where i?io = i?/10^°cm (see e.g. Bchreibcr & Gansicke 
2001|; If ameury et al. 1998|; [Ludwig et al. 1994D p|. 



Taking into account heating of the outer edge of the 
disc by the impact of the stream, the critical mass transfer 
rate is given by an implicit equation as the turning point 
of the modified S-curve Ma,! depends on the mass trans- 
fer rate itself; Ma,! = MA,i{Rio, Afwd, Mtr)- The critical 
mass transfer rate is then given by a solution of 

Mt, = MaARio, Mwd, Mtr), (7) 

(see also ^tehle et al. 2001 ) . Note that because the stream 
impact heating decreases with decreasing radius, the crit- 
ical mass transfer rate is not necessarily defined by the 
vertical structure at the outer edge of the disc. In order to 
solve Eq. (7), we modify the approximation of the upper 
turning point of the S-curve by including impact heating. 
We obtain for the parameter range of SS Cyg: 



log(MA,i) = log(MA) - 0.039 



Qi 



IQiOergs-i' 



(8) 



Using this description, we solve Eq.(6) and obtain a series 
of solutions. The results can be approximated with the 
following power law: 
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Merit,i = 3.18x 10" gs 



R 



10 



5.83 



2.750 



vi.oy 



-0.267 



1.19 Afp 



-0.287 



This leads to 



(3.2 ± 0.6) X IG^^gg 



(9) 



(10) 



Fig. 4 shows S-curves for the outer edge of the accre- 
tion disc assuming the above obtained mean mass trans- 
fer rate and considering the additional heating by the 
stream (Eq. (2)). The transition between ahot — 0.1 and 
Q^coid — 0.02 is assumed to be given by the following tem- 
perature dependence of a: 



log(a) = log(Q;coid) + [(log(Q!hot) - log(Q;coid 

'2.5 X lO^K 



(11) 



(see Hameury et al. 1998). 

For the mean mass transfer rate we obtain Mtr ~ 2 x 
Ma,! (Fig. 4), which would lead to steady disc accretion. 
Considering the uncertainties of the involved parameters, 

we find Mtr^Mcnt,i- 

4. The onset of the cooling front 

The critical and the mean mass transfer rate derived in the 
previous section are sensitive to the uncertainties in the 



The wea k dep endence of Ma on Qh obtain ed by Hameury 
et al. (1998 ) and Schreiber fc Gansicke (2001 ) is negligible in 
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Figure 4. S-curves for outer radii [in units of 10^° cm, 
R = 0.9i?Li = 5.83 (dotted), 0.8 Rli =_5.18 (solid), 
0.7 Rli = 4.53 (short dashed hne)] using Mtr = 4.2 x 
lO^^gs"^ (solid horizontal line, see Eq.(5)). At the outer 
edge of the disc (dotted line) the instability is completely 
suppressed by the additional heating due to the stream 
{iji = 1.0). Moving inward, the additional heating de- 
creases (see Eq. (2)) but the mean mass transfer rate (solid 
horizontal line) is still clearly above the value critical for 
stability. The critical accretion rate Afcrit.i (see Eq. (10)) 
is represented by the dashed horizontal line. Notice, the 
S-curves are calculated with Mti — Mtr 7^ Mdit.i- 

binary parameters and the assumed efficiency of stream 
impact heating rji as the additional heating can signifi- 
cantly reduce the critical mass transfer rate. Thus, while 
our results presented in Sect. 3 suggest that the DIM fails 
in describing the observations of SS Cyg, they do not allow 
to categorically exclude the validity of the DIM. 

However, in this section we show in an independent 
way that the observed light curve of SS Cyg provides ad- 
ditional evidence for the failure of the DIM, independent 
from the details of stream impact heating or the uncer- 
tainties in the system parameters. For the mass accretion 
rate at the onset of the decline Afotd the DIM predicts an 
essentially lower value than required by the observations 
as, (a) the DIM predicts that the start of the cooling front 
at the outer edge of the disc defines the onset of the decline 
in the outburst light curve. The cooling front starts when 
the accretion rate of the quasi stationary disc is exactly 
equal to the critical rate, thus. 



M, 



otd.DIM 



Mr: 



(3.2 ± 0.6) X 10 



17. 



(12) 



(b) the light curve of SS Cyg shows that there is no no- 
ticeable difference between the absolute magnitude at the 
onset of the decline and the outburst magnitude rriy = 
8.5 ± 0.1 (Fig. 1) and therefore, the observed light curve 
constrains the accretion rate at this point to be equal to 
the outburst accretion rate 



the context of this paper. 



M, 



otd,obs 



Mo 



(3.62 ± 1.47) X 10^** gs 



(13) 
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The obtained huge discrepancy between the prediction of 
the DIM for Motd (Eq. (12)) and the value required to 
explain the observed visual magnitude (Eq. (13)) clearly 
confirms the failure of the DIM. 

It is worth noting that even if we completely neglect 
stream impact heating the DIM cannot explain the obser- 
vations. Then, the critical mass transfer rate can be ap- 
proximated by Eq. (6) and the accretion rate necessary to 
reach the visual flux during outburst somewhat increases: 

Motd^DiM = Afcrit = (7.6 ± 2.0) X 10i^gs-\ (14) 
Motd^obs - Mout = (3.71 ± 1.51) X 10i"gs-\ (15) 
Mtr = (4.3 ± 1.8) X lO^^gs^^ (16) 

Obviously, the problem for the DIM obtained in Sect. 3 
vanishes if we neglect stream impact heating of the outer 
disc. Nevertheless, the DIM clearly fails as Motd.obs ^ 

-Motd, DIM- 

5. Discussion and Conclusion 

Considering the "larger" distance of SS Cyg and the ob- 
served long term light curve, we calculated the mean mass 
transfer rate Mtr- Taking into account stream impact 

heating of the outer disc we find Mtr ~ 2 x Merit (Fig- 4) 
which contrasts with the predictions of the DIM: with such 
a high mean mass transfer rate, SS Cyg should not be a 
dwarf nova at all, but a novalike variable with a stationary 
hot accretion disc. This result is, however, sensitive to the 
assumed efficiencies of stream heating and disk irradiation 
and the uncertainties of the binary parameters, the incli- 
nation, and the distance. While this result on its own is 
not sufficient to proclaim a general failure of the DIM, an 
independent suggestion for a major problem with the DIM 
comes from comparing the accretion rate predicted by the 
DIM at the onset of the decline with the accretion rate 
required by the observed visual magnitude at this point 
(see Sect. 4). Hence, the conclusion from both these results 
must be; the current DIM cannot explain the observations 
of SS Cyg for a distance of 166 ± 12 pc. 

There are three possibilities that may resolve this prob- 
lem: (a) the temperature dependence of a, the key ingre- 
dient in the accretion disc limit cycle model, has to be 
modified, (b) increased mass transfer during the outburst 
plays an important role, or finally (c) the HST/FGS dis- 
tance of 166 pc and/or the assumed inclination are wrong. 

Considering (a), it seems that we do not yet have 
an adequate understanding of the viscosity in accretion 
discs. The instability necessary to obtain the limit cycle 
behaviour appears mainly due to an artificial change of 
a, given by Eq. (11), at the point where the small "natu- 
ral" , i.e. without a change of a, ionisation instability has 
been found in the S-curves. There is no physical reason for 
this assumption and it might be the case that a starts to 
change at higher temperatures than previously assumed. 
This would lead to a higher value of the critical mass trans- 
fer rate. As such a modification would affect the dividing 



line between dwarf novae and nova-likes, one should care- 
fully consider the consequences for other dwarf nova before 
changing Eq. (11). 

With respect to (b) , we point out that another possibil- 
ity to solve the problem obtained in Sect. 2 - 4 is to assume 
that the mass transfer increases during the outbursts, e.g. 
due to irradiation of the secondary star. Assuming a mean 
mass transfer rate of Mtr, qui = 1-5 x 10^''gs~^ during 
quiescence, rise, and decay, the mean mass transfer rate 
during the outburst state should be 



Mt 



Mout tqs — {tqui + ^ris + tdcc)M 



tr,qui 



^out ^ris ^de 



= 2.4 X 10^^gs"\ 



(17) 



i.e. increased over the mean quiescent value by a factor of 
> 15. In this scenario the onset of the decline in the light 
curve would represent decreasing mass transfer and not 
the start of the cooling front. 

Finally, considering (c), we can only note that at least 
one of the additional CV parallaxes measured with HST, 
namely U Gem, agrees very well with previous distance 
estimates (Harrison et al. 1999). UGem has a notice- 
ably lower accretion rate than SS Cyg, therefore the dis- 
tance estimated based on K band magnitude and/or spec- 
troscopy of the secondary star is much less affected by 
emission from the disc than in SS Cyg. Thus, at least 
for UGem an error in the HST/FGS parallax appears 
extremely unlikely. However, assuming i = 37° and de- 



manding A/o 



1.0 X IQi^gs-i for SSCyg leads to 



Aftr ~ 1.2 X lO^^gs ^ and requires (i ~ 1 17 pc for agree- 
ment with the DIM. 
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